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Abstract — A limited network of four dual-channel microwave 
radiometers^ with frequencies of 20.6 and 31 65 GHz, was operated 
in the front range of eastern Colorado from 1985 to 1988. ua 

from November 1987 through October 1988 f ® an ^ Y ^°J^e L 
both single-station and joint-station brightness temperature and 
attenuation statistics. Only zenith observations were ™de- The 
snatial separations of the stations varied from 50 km to 19U Km. 

Before ^statistics were developed, the . dat ^f e t ^ e ®f C 2o1 
rigorous quality control methods. One such method, that of 20.6 
vs. 31.65 GHz scatter plots, is analyzed in detail, an 
comparisons are made of measured vs calculated 

65 GHz vertical attenuations of 5 and 8 aB are exceeae 
of 'the time. For these four stations and at the same (K01% level 
diversity gains from 6 to 8 dB are possible with the 50 to 190 km 

separations . 

I . introduction 

During 1985 through 1988, the Wave Propagation Laboratory 
/ W pl) of the National Oceanic and Atmospheric Administration 
NO^i) operated a limited network of four dual-ch 
radiometers. The radiometers, operating at 20.6 and 31.65 GHz, 
were deployed for the meteorological purposes of measuring 
orecipitable water vapor (PWV) and integrated cloud liquid water 
?CLW)^ Data from ground-based microwave radiometers are commonly 
used to derive attenuation statistics (Ortgies et al., 1990, 

Vogel etal! 1991; Fionda et al., 1991). Using the radiometric 
technique, these authors have shown that attenuation levels JP 
about 12 dB can be derived with good accuracy. Since attenuation 

?sf he sum of absorption and scattering and rf ?fmf hr ?eg?!n 
not negligible when rain rates are in the 5 to 10 mm/hr region, 
rad iometr ica 1 ly-der i ved absorption is not representative of 
attenuation under such rain conditions. Data from two of the 
Colorado stations, Denver and Platteville, have previous y 
used to derive single-station and two-station attenuation 
s?a?is?ics fo? 2 three-month periods during 1988 ( Fionda et al., 
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1991) . Here, we greatly extend that work by deriving attenuation 
and attenuation-diversity statistics for an entire year's data 
taken at each of the four stations. For completeness, we also 
present complementary statistics on the basic measured variable - 
brightness temperature T b . Since operations at two of the 
stations, Fleming and Flagler, were discontinued during November 
of 1988, we used data starting from November 1, 1987, and ending 
at October 31, 1988, to derive a complete year's statistics. 

II. Description of Radiometers in the Colorado Research Network 

A. Location 


A description of Wave Propagation Laboratory radiometers was 
given by Westwater and Snider (1987); here we will only describe 
in limited detail the dual-channel radiometers of the Colorado 
Research Network. Their geographical coordinates and the 
separation distances are given in Table 1. 

B. Radiometer Characteristics 

The dual-channel radiometers were designed, constructed, and 
field-tested by WPL; a complete description of the systems is 
given by Hogg et al. ( 1983) . The instruments were designed to run 
continuously, to provide unattended operations, and to operate in 
almost all weather conditions. The salient characteristics of 
the instruments are shown in Table 2. Note that the antenna 
beamwidth at Denver differs from those of the other three 
stations. Field experiments (Snider, 1988), in which a steerable 
radiometer with a 2.5 deg beamwidth was compared with the network 
radiometers with their 5 deg beamwidths, showed a 0.99 
correlation between the systems. The receivers of all four 
radiometers are of the same construction. The internal 
calibration of the radiometers is done by switching between the 
antenna and two hot blackbody loads; external calibration is done 
approximately every two weeks using the "tip cal" method (Hogg et 
al. , 1983; Decker and Schroeder, 1991). 

c * Methodology to Derive Attenuation from Emission 

The basic quantity measured by a radiometer is brightness 
temperature, which is closely related to input power present at 
an antenna (Ulaby et al., 1981). Although the probability 
distributions of T b are of interest in themselves, the quantities 
needed by communication engineers are the distributions of 
attenuation. We derive attenuation t (dB) from T b by using the 
well-known formula (Westwater et al., 1990) 

T (dB) = 4.34 ln{ (T m - T c ) / (T m - T b ) } , ( 1 ) 

where T m = mean radiating temperature (K) , 

and T c = cosmic background temperature = 2.75 K. 
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In deriving t, we used monthly mean values of T. (see Table 3) 
that were calculated from our radiative transfer and cloud 

models. With monthly mean values of T„, est i m ^ e ^ h ^ S d ®^°in q 
about 7°C occur. These errors become important when derivi 9 
at?enu!tion from the higher values of T b , say those greater than 

150 K. 

III. Quality Control 

The radiometric data were taken by radiometers ?P®* ated 

in an unattended mode, although bimonthly on-site ? a ^”f^ onS 
were done. For the most part, the data were of high quality, 
although occasional outliers had to be removed from the data. 

Such outliers can arise from liquid and ice buildup on . 

antennas spurious signals of electromagnetic origin, calibration 
-f t-c in the receivers, and data transmission errors. To 
SSSUS obvious 6 erroneous data, we plotted and inspected daily 
time series of the following quantities: brightness temperature 
i b at 20.6 and 31.65 GHz; derived values of PWV and CLW . If a 
record had an obvious error at either frequencies.data fr th 
entire record were removed. Next, scatter plots of T b s at born 
frequencies were constructed; usually, suspicious points were 
easily identified from these plots. Westwater and Falls (1991) 
described the method in more detail and gave examples of the 
method . That method was applied to all the data analyzed m this 

report . 

One other consideration is necessary before we P re sent our 
results Westwater et al. (1990) have shown that calculations of 
clear air brightness temperature, based on rad *° S ° nde 

soundings and contemporary absorption algorithms (hieb < e a n te d 

t avton * 1987 } differed from T h measurements that were calioratea 
by V the'tip cii method. Consequently, we _ adjusted our calculated 


by tne tip cai iuu. , --- 

values T^ ,, based on their absorption algorithms, to be 
consistent with those determined by a tip cal procedure, 
adjustedvalues T b adJ , based on the data set described by 
Westwater et al. (1990), are 


These 


20.6 GHz: 
31.65 GHz: 


Tb.adj = 1-144 T b>cal - 0.049 
T b .adj = °- 970 T b,cal + 1 - i01 


( 2 ) 


1 . 


A scatter plot of adjusted calculated data is shown in Fig 
Se nSte that the range of calculated T b ’s is considerably less 
i-han that of the measured values, presumably because we 
underes?U£e £he range of cloud liquid that occurs This is not 
surorising, because our cloud calculations, based on the Rayleig 
aDoroximat ion , do not contain a rain model. However, over the 
- i _..-i t •« the slopes and ranges of variation of 

measured SS°£E3a&- are "similar. To gain insight into 

the ranges of variation and their causes, we also calculated T b 
as a function of (a) PWV for clear conditions, (b) CLW tor cio y 
conditions, (c) PWV during both clear and cloudy conditions, 

ClS for both clear and cloudy conditions. These results are 
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shown in Figs. 2 and 3 for 20.6 and 31.65 GHz. We note from 
these figures that for clear conditions, T h 's at both frequencies 
vary linearly with PWV over the range of o'. 3 to 4.0 cm. The 
interesting behavior occurs when clouds are present; for a given 
vapor, there is a considerable range of variation in T K due to 
JJ® allowable range of CLW. Conversely, for a given amount of 
clw, there is a considerable range of variation in T h due to 
clouds. These figures also clearly show the relative sensitivity 
of the two channels to vapor (20.6 GHz is more sensitive) and to 
liquid (31.65 GHz is more sensitive). Finally, we show in Fig. 

4/ ®?fV ter Plots of calculated T b 's, during clear and cloudy 
conditions. Thus, the lower straight line borders in the scatter 
plots are determined by the amount of PWV, and the more diffuse 
upper boundaries are determined by the amount of CLW. 

Iv * Single- Station Attenuation Statistics 

After the quality control methods discussed above were 
applied to the data, cumulative distributions of brightness 
temperature and radiometrically-derived attenuation were derived 
for each of the four stations. In addition, the composite of all 
stations was also computed. The results for brightness 
temperature and attenuation are shown in Figs. 5 and 6. Perhaps 
not too surprisingly, there is not a substantial difference 
between any of the four stations. At both frequencies, Denver is 
slightly colder in T b and is slightly less attenuating than the 
other three stations. For both frequencies, attenuation is less 
^ h ^„ 1 . dB at the 1% level, and is less than 7.5 and 9.5 dB at the 
0.01% level for 20.6 and 31.65 GHz. Finally, the four-station 
composite values are plotted in Fig. 7. It will be of interest 
to compare these radiometrically derived values of attenuation 
with forthcoming beacon measurements from ACTS (Chakraborty and 
Davanan, 1991) . 1 

V* Joint- Station Attenuation Statistics 


At the time these data were taken, the sampling times for 
each station were not synchronous, and starting times could 
differ by up to 1 1/2 minute. To compute joint-station diversity 
statistics, it is necessary to put time series from each station 
into one-to-one temporal correspondence. For our data, we set up 
a 2-minute window, and when starting times from each of the two 
stations fell within this window, the data were placed in 
correspondence. If a complete pair of data was not available 
w? 2 “ minute sample was eliminated. Over a year's time, we were 
3 s f9 ni 5 i( f ant sample size for analysis: a minimum 
~ ? 39 data . pairs (Platteville - Flagler) and a maximum of 

218,521 data pairs (Denver - Fleming). 


. *^ e results for the joint-station diversity analyses are 

given m Fig. 8. We note that the diversity curves are all quite 
similar for the stations that are separated by » 150 km, but that 
the closer pair, Denver - Platteville, differ significantly from 
the other five station pairs. Roughly, at the 0.01% level, the 
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• c - n on „er - piatteville are about 6 dB at both 
diversity gains for Denv 150k i separations, the corresponding 

gains^re aiout "-7 and 5 dB at the lower and upper frequency. 

VI. summary and Plans 

We have derived yearly cumulative ^^“foulTs^tions of 
brightness temperature and at tenuati. “iJigle aid joint-station 
the Colorado Research Network. stations whose separations 

distributions were derived for these stations det £ iled 

iaiarsu--. 

Ssir^U^i at°Norman , Oklahoma , and will soon be processing 
these data. 
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Table 1. 


Geographical coordinates and spacings of the Colorado 
Research Network 

Site Latitude (deg. ) Longitude (deg. ) Elevation (msl ) 


Denver 

39.75 

N 

104.87 

W 

Flagler 

39.12 

N 

103.09 

W 

Fleming 

40.63 

N 

102.94 

W 

Platteville 

40.18 

N 

104.73 

W 


1.611 

1.463 

1.337 

1.523 


Combination 


Spacing (km) 


Denver-Platteville 49.25 
Denver-Flagler 168.04 
Denver-Fleming 190.79 
Flagler-Fleming 168.28 
Flagler-Platteville 183.19 
Fleming-Platteville 159.50 


Table 2. Characteristics of WPL Dual-channel Radiometers 


Operating frequencies (GHz) 
Viewing 

Antenna half-power beam width 

Total bandwidth 
(double side band) 

Integration time 

Sensitivity (for 2-min 
integration time) 

Estimated absolute 
accuracy 


20.6 and 31.65 GHz 
Zenith 

2.5° (Denver), 5.0° (other 
three) 

1 GHz 

2 min 

0.05 K rms 
0.75 K 


Table 3. Monthly-averaged mean radiating temperatures/standard deviations 
for 20.6 and 3L66 GHz., Denver, Colorado, 1970 - 1986. 















Month 

1 

2 

3 

4 

6 

6 

7 

8 

9 

10 

11 

12 


20.60 269.7 260.3 262.8 266.4 270.7 276.4 278.3 277.7 274.3 270.9 264.4 260.8 

6.6 6.3 6.4 6.7 6.2 6.4 8.9 4.3 6.9 7.4 6.6 6.8 

8L66 266.7 266.3 268.6 262.1 266.3 270.9 278.8 273.3 269.8 266.4 260.1 266.7 

6.0 6.0 6.4 6.7 6.1 6.6 4.4 4.7 6.9 7.2 6.0 6.2 
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AdJ. Tb at 20.60 GHz (K) 

1 * s ^ atter Plots of measured (A) and calculated (B) T at 
20.6 and 31.65 GHz from Denver, Colorado . The measured L, 
were taken in 1988; the calculated data were based on rldioJondes 
taken during 197o - i9 85 . The calculated T b values have been 
adjusted to be consistent with those determined by the tip cal 
procedure (see Section III for details) . P 
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Fig . 2 . Scatter plots of ad j u ste d ca 1 culated T at 2 0 * 6 
fa\ PWV for clear conditions; (b) CLW tor ciouuy i 
c p!!v fir both clear and cloudy conditions,; and (d) CLW for 
both clear and cloudy conditions. 
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20.60 (GHz) Zenith Allen. (dB) 20.(0 (GHz) Zenith Bright. Temp. (K) 



*-• Oil* s 1* *0 30 40 00 00 70 00 10 0 } 10 00 00.0 00.0 00.00 


Percentage of Time that Ordinate is Exceeded 

Fig. 5. Cumulative brightness temperature (A) and zenith 
attenuation (B) statistics for the four stations of the Colorado 
Research Network at 20.6 GHz. 
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31.65 (GHz) Zenith Brijhl. Temp. (K) 


11 / 1*87 - 10 / 1*88 


-V 

A\ Denier 

\\\ PltUeville 

Vi ’ Flemioi 
V ’ FU|Ur 

V\ 'Total 


Mean Tb(K) Std. DeT.(K) Ob». 


— 14.033 

- - 15.534 

11.078 
-- 18.358 
— 15.472 


*.081 

10.213 

12 

*.*46 

10.444 


23 ) 52 * 

237740 

245878 

204701 

*27848 


\ 


, * . .8 I. I. «. .0 •• T. .0 .« •» •••••• •»• • 


11 / 1*87 - 10 / 1*88 


Deafer 

Plitlefille 

fleniai 

FUfUr 

Total 


Mean Aban(dB) Std. Dev.(dB) Obs. 


0 . 1*5 

0.22 

. . . 0.231 

0.234 

0.22 


239518 

237724 

245810 

2048*3 

927743 


# I I III II I I. I I l ■ ‘ ' 

... , . it " 

Percentage of Time that Ordinate is Exceeded 

su.irr sasars 

Research Network at 31.65 GHz. 




Zenith Alien. (48) Zenith Bright. (K) 


300 



a F t?enu;tio“(B P ritatistSs n ?or ^2 P ? rature < A > «nith 
Research Network. th four statlons of the Colorado 
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Zenith Atten. (dB) 


Max / Min ( Denver .Platteville) 


t 1 / 1917 - 10/IISI 


Main Aban(dB) Sid. D«T.(dB) Oba. 


__ Ui 1.532 

20.60 ui* . - - 0.211 

r 3i.65 i" — -- 

\ Mia 0.181 


0.274 110145 
0.139 101145 
1.374 111133 
9.111 110033 







Zenith Allen. (dB) Zenilh Atten. (dB) Zenith Atten. (dB) 
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